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ABSTRACT 

We use 62,185 quasars from the Sloan Digital Sky Survey DR5 sample and standard 
virial mass scaling laws based on the widths of H/3, Mgll, and ClV lines and adjacent 
continuum luminosities to explore the maximum mass of quasars as a function of 
redshift, which we find to be sharp and evolving. This evolution is in the sense that 
high-mass black holes cease their luminous accretion at higher redshift than lower-mass 
black holes. Further, turnoff for quasars at any given mass is more highly synchronized 
than would be expected given the dynamics of their host galaxies. We investigate 
potential signatures of the quasar turnoff mechanism, including a dearth of high-mass 
quasars at low Eddington ratio. These new results allow a closer examination of several 
common assumptions used in modeling quasar accretion and turnoff. 

Key words: black hole physics — galaxies: evolution — galaxies: nuclei — quasars: 
general — accretion, accretion discs 



1 INTRODUCTION 

Supermassive black holes (SMBH) are found today at 
the centers of nearly every galaxy where there have 
been sensitive searches. These SMBH are the remnants 
of bright quasars which are found in many galaxies at 
redshifts z ~ 2. The mechanisms involved in the quasar 
turnoff process are not well understood (cf. Thacker et al. 
2006). In this work, we combine the large quasar sam- 
ple provided by the Sloan Digital Sky Survey (SDSS) 
([Schneider et al. 2007[) with recent advances in virial mass 
estimation (Vestergaard & P eterson 2006|IShen et al. 2008|) 
to develop new constraints on quasar turnoff. 

It had been previously thought that most quasars were 
within a factor of ~ 10 of their Eddington luminosity 
(Kollmcier ct al. 2006), and quasars have been commonly 
modeled as 'light- bulbs', operating either near Eddington 
or without substantial luminous accretion, although this 
is likely an oversimplification ([Hopkins 2008]). In Paper I 
USteinhardt fc Elvis 2 009 ) , we used the quasar distribution 
in the mass-luminosity plane as a function of redshift. Had 
the relationship between quasar mass and luminosity been 
trivial (as implied by the the light-bulb approximation) the 
mass-luminosity plane would not present an improvement 
over the quasar mass function or luminosity function. In- 
stead, we found a more complex relationship between mass 
and luminosity than had been previously assumed. Several 
new puzzles have emerged from this view of the quasar dis- 
tribution, including a sub- Eddington boundary that restricts 
luminosities of the highest-mass quasars at every redshift to 



lie below their Eddington limit. We also noted an upper 
mass limit varying with redshift. In this work, we use this 
new, multi-dimensional view of the quasar distribution to 
investigate this high-mass 'turnoff' in detail. 

We use these results to evaluate three basic assump- 
tions: (1) that the dynamics of quasar turnoff are closely 
linked to the dynamics of the host galaxy, (2) that a seg- 
regation of quasars by mass is equivalent to a segregation 
of quasars by luminosity, (a corollary of the light-bulb ap- 
proximation), and (3) that SMBH undergo no fundamental 
change during turnoff, but rather run out of fuel (which im- 
plies that an individual SMBH might have many luminous 
episodes during its lifetime and would undergo another one 
today given sufficient fuel). 

In § [2l we briefly review the SDSS sample and exhibit 
its behavior in the M — L plane. A full discussion of this 
sample and its uncertainties can be found in Schneider et 
al. (2003) and Paper I. In § [3] we establish that massive 
quasars are indeed turning off, and doing so earlier than 
lower-mass quasars. In § 21 we use this to select individual 
quasars near turnoff and search for signatures of their turnoff 
mechanisms. In § O we show that the synchronization of 
turnoff of different quasars at a given mass is closer than that 
of the dynamics of their host galaxies. We also show that the 
synchronization of quasars at fixed mass is characteristically 
different from that of quasars at fixed luminosity as reported 
in Amarie et al. (2009). In § [6l we evaluate the three basic 
assumptions above in light of this new evidence. 
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Table 1. Summary statistics on quasars in the 13 redshift and 
emission line bins. 
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Figure 2. The SDSS quasar locus of the Hj3 mass sample in 
the M — L plane at redshift 0.2 < z < 0.4. The locus should 
be bounded by SDSS detection limits, Lfidd (dashed line), 
and on the high-mass end by an unknown mechanism responsi- 
ble for quasar turnoff. There is an additional sub-Eddington 
boundary (Steinhardt & Elvis 2009) with slope below that of 
^Edd- Finally, high-mass, low-luminosity quasars are sparse, sug- 
gesting the possible presence of an additional boundary (HM- 
LLB) or added complexity in the nature of quasar turnoff. The 
bright-object SDSS saturation limit does not intersect the 
quasar locus. These lines are defined quantitatively either here 
or in Paper I. 



3 MASSIVE QUASARS ARE TURNING OFF 

As shown in Table [1] both the maximum luminosity and 
the maximum mass of the quasar locus are monotonically 
decreasing towards lower redshift. Figure [T] hints the same 
might be true of the average luminosities and masses of 
the quasar populations in each bin. While it has been 
observed that individual quasars have variable luminosity 
| |Mathews~fc Sandage 1963]), we know of no mechanism by 
which the central black hole can substantially reduce its 
mass. Therefore, we must interpret Figure [T] as showing us 
that the most massive quasars in any cosmological epoch are 
in the midst of disappearing. 

Each redshift bin in[T]is of an equal size in redshift, and 
therefore a different size in comoving volume. It is therefore 
expected that the most massive quasars become rarer at 
low redshift; Mgll-based mass estimates at 1.8 < z < 2.0 
encompass a comoving volume 15 times larger than Mgll- 
based mass estimates at 0.6 < z < 0.8. If the number density 
of massive quasars were constant, there should be 15 times 
as many high- mass quasars visible at 1.8 < z < 2.0 than at 
0.6 < z < 0.8. 

The number density at a given mass as a function of 
redshift for all quasars in the SDSS catalog is shown in Fig- 
ure [3] In all but the lowest mass bins, the catalog comoving 
number density rises to a peak between a redshift of 1 and 
2, then declines towards lower redshift. However, the lumi- 
nosity distribution for quasars at fixed mass and redshift is 
cut off on the low-L end by the SDSS detection limit for 
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2 SDSS QUASARS IN THE M-L PLANE 

Black hole masses for 62,185 of the 77,429 SDSS DR5 
quasars (Sch neider et al. 2007P were determined by Shen 
et al. (2008) using H/3- and ClV-based virial mass esti- 
mates from Vestergaard & Peterson (2006) and Mgll-based 
estimators from McLure & Dunlop (2004). As discussed 
elsewhere l|Shen et al. 2008] |Risaliti, Young, fc Elvis 2009| 
ISteinhardt fc Elvis 20091) , ClV-based estimation is substan- 
tially less accurate than H/3 or Mgll, while Mgll has a 
smaller statistical uncertainty than H/3. Vestergaard & Pe- 
terson (2006) report an uncertainty for virial mass estima- 
tion of ~ 0.4 dex from calibration against reverberation 
mapping, while in Paper III we discuss the possibility that 
the statistical uncertainty may be closer to ~ 0.15 dex and 
the remaining deviations due to systematic effects. 

As in Paper I, we divided the sample into redshift bins 
of size 0.2 in z for z < 2.0 and larger redshift bins for z > 2.0. 
Table [T] contains summary statistics on the objects in each 
bin. Within each bin, there is a distribution of black hole 
masses spanning ~ 2 dex (Figure [TJ. 

Figure [T] shows quasar loci in the L — M plane for each 
of the 12 redshift bins as contour plots. An individual quasar 
locus in Figure [2] has been labeled at its various boundaries. 
The quasar locus at each redshift is bounded by SDSS detec- 
tion limitations (bottom) , the Eddington luminosity (dashed 
line, upper- left), a sub-Eddington boundary (top) discussed 
at length in Steinhardt & Elvis (2009), a high-mass cut- 
off (right) associated with quasar turnoff and discussed in 
§ and possibly an additional high-mass, low-luminosity 
boundary discussed in § 2] 

No quasars in the catalog have luminosity greater than 
LEdd with statistical significance, while a sub-Eddington 
boundary is present in each panel. In each panel, the low- 
luminosity boundary on the quasar locus is the only bound- 
ary due to SDSS detection limitations. In particular, SDSS 
would be capable of detecting quasars past the high-mass 
and high-luminosity boundaries of the quasar locus in each 
panel and including them in the catalog if they were to exist. 
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Figure 1. A contour plot of the black hole mass-bolometric luminosity distribution in 12 different redshift ranges, as described in 
Steinhardt & Elvis (2009). In each redshift bin, the quasar number density has been normalized to the peak number density. The top 
three panels are from the low-redshift sample, i.e., with H/3-based mass estimates; the middle six panels are from the medium-redshift 
sample, i.e., with Mgll-based mass estimates; and the bottom three panels are from the high-redshift sample, i.e., with CIV-based mass 
estimates. The red dashed line is drawn at the sub-Eddington boundary (Paper I) and the blue dashed line in the Mgll panels are drawn 
at the high-mass, low-luminosity boundary. The green dashed line is drawn at approximately 19.1 in i band for the average redshift in 
each panel; below this line objects are from the serendipitous sample. 



many masses (e.g., at 9.0 < log M/M Q < 9.25 in Figure g]). 
Therefore, the total number density of quasars at a given 
mass is only known in some redshift bins, for which only 
a negligible fraction of quasars at that mass lie below the 
SDSS limit. To address this bias, we only include bins for 
which the number density at the detection limit is < 20% 
of the peak. Figure U shows a sample cut for quasars with 



9.0 < log M/Mq < 9.25. We include several bins where there 
are fewer than ten total quasars, but none at the detection 
threshold. Restricting Figure|3]to those points produces the 
number density-lookback time plot in Figure [S] 

The number density distributions in Figure[S]clearly de- 
cline with increasing cosmic time (decreasing redshift). The 
higher-mass curve appears to decline faster. To investigate 
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Figure 3. Raw comoving number density evolution of quasars 
in four different mass bins at 0.2 < z < 3.2, normalized to 
peak. Poisson uncertainties are indicated and the comoving num- 
ber density at each mass is normalized to its peak value. Each 
curve contains quasars in a 0.25 dex range in mass: green 9.0 < 
log M/M e < 9.25, cyan 9.25-9.5, blue 9.5-9.75, and magenta 9.75- 
10.0. 
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Figure 4. Luminosity distribution of SDSS quasars with masses 
9.0 < log M/Mq < 9.25 in different redshift bins: red 0.8 < z < 
1.0, yellow 1.2 < z < 1.4, green 1.6 < z < 1.8, cyan 2.0 < z < 2.2, 
and blue 3.0 < z < 3.2. Poisson uncertainties are indicated and at 
each redshift the number density is normalized to its peak value. 
The red, yellow, and green curves dip below 20% of peak at the 
detection limit and therefore are considered sufficiently complete, 
while the cyan and blue curves are not. 
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Figure 5. Completeness-selected comoving number density evo- 
lution with lookback time for quasars in different mass bins at 
0.2 < z < 2.0. Displayed data is restricted to bins in which 
nearly all quasars could be detected and included in the SDSS 
catalog. Poisson uncertainties are indicated and the comoving 
number density at each mass is normalized to its peak value. 
Each curve contains quasars in a 0.25 dex range in mass: green 
9.0 < log M/Mq < 9.25, cyan 9.25-9.5, blue 9.5-9.75, and ma- 
genta 9.75-10.0. 



Table 2. Best-fitting exponential decays for declining comoving 
number densities of quasars in different mass bins. Also included 
are the best-fitting declines using only Mgll-based mass estima- 
tion. 
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this possibility, we fit this decline as an exponential decay at 
each mass. The results are summarized in Tableland Fig- 
ure El At higher masses, characteristic time of the decline 
in comoving number density is indeed shorter. The higher 
mass uncertainty in the H/3-based mass estimates leads to 
a better fit when only the Mgll-based mass estimates are 
used, but both estimates are in agreement at each mass. The 
high x 2 /DOF values either suggest that the decay cannot be 
well-modeled as exponential or that Poisson errors are un- 
derestimates of the true uncertainty in the comoving number 
density. The latter almost certainly applies, as the mass un- 
certainty is ~ 0.4 dex (Vcstergaard & Peterson 2006) while 
bins are just 0.25 dex in width (although in Paper III we dis- 
cuss evidence that the statistical uncertainty in virial mass 
estimation may be smaller than previously believed). How- 
ever, the decline may not be exponential. An improvement 
in the detection limit would allow a measurement of the co- 
moving number densities over a wider range of redshift, and 
therefore might allow a better determination of the exact 
nature of their declines. 
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Figure 6. Timescale for best-fitting exponential decay in comov- 
ing number density of quasars as a function of mass. The black 
points have been fit using both H/3 and Mgfl masses, while the 
red points only use Mgfl mass estimates. 



The timescales for decline Td are consistent with a drop- 
ping exponential with mass (Figure [6}. Using all mass es- 
timates available, r d oc M" 2 ' 097 * ' 004 with a xV D OF of 
0.60, while the timescales using only Mgfl mass estimates 
are proportional to r d oc M" 1 ' 809 * 155 with a x 2 / D °F of 
0.626. The fit uncertainty using all mass estimates is small, 
but the uncertainties in individual measurements of Td indi- 
cate that the true uncertainty in slope is considerably higher 
than the fit uncertainty. 

The comoving number density of high-mass quasars 
(Figure[5]) shows that high-mass quasars are indeed less com- 
mon at lower redshift, with the highest-mass quasars having 
the most rapid decline in number density. At each redshift, 
the SDSS catalog contains a snapshot of a population of the 
highest-mass quasars set to decline. Massive quasars indeed 
turn off, with the highest-mass quasars turning off first. In 
§ [4] we search for signatures of this quasar "downsizing' in 
an effort to study the causes and mechanisms involved. 



4 SIGNATURES OF QUASAR TURNOFF 

At every redshift, the high-mass quasar population as a 
whole is declining rapidly as in § [3] It follows that a snap- 
shot of that population must catch many individual high- 
mass quasars that will shortly undergo turnoff. Further, if 
the turnoff process is gradual, it is very likely that the SDSS 
catalog will contain quasars in the midst of turnoff. For ex- 
ample, if quasar turnoff takes a minimum of 40 Myr as sug- 
gested by Hopkins et al. (2008), we should expect ~ 200 
quasars with Mbh > 9.5 at 1.8 < z < 2.0 to turnoff in the 
next 1 Gyr, and therefore > 8 to be in the midst of turnoff 
while emitting the light we observe. 

Here, we search for anomalous properties of the highest- 
mass quasars in each redshift bin the hopes that they will 
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Figure 7. The high-mass, low-luminosity end of the quasar mass- 
luminosity plane in four redshift bins. Quasars at 0.4 < z < 0.6 
are displayed with H/3 mass estimates, while quasars at 0.8 < 
z < 1.0, 1.2 < z < 1.4, and 1.6 < z < 1.8 use Mgll mass esti- 
mates. Quasars brighter than i = 19.2 are shown in green, while 
quasars in the serendipitous sample are shown in red. High-mass, 
low-luminosity quasars are scarce in three panels, indicating the 
presence of an HMLLB (blue). 

provide signatures of the mechanism by which quasars turn 
off. 



4.1 The High-Mass, Low-Luminosity Boundary 

We can attempt to probe the turnoff rate by looking at the 
luminosities of the most massive quasars. Consider one of the 
most massive quasars in a given redshift bin that lies at a lu- 
minosity near the sub-Eddington boundary (SEB, Paper I). 
As discussed above, this is a snapshot of a quasar which will 
shortly begin its turnoff. If this turnoff is gradual, a snap- 
shot of this quasar at a later time must show a mass larger 
than any quasar at the SEB but at a luminosity well below 
the SEB. Since Mbh for the most massive quasars declines 
with decreasing redshift, a more gradual turnoff process will 
result in a larger population of massive, low-Eddington ratio 
quasars with mass more greatly exceeding the most massive 
quasars lying on the SEB. 

However, Figure [T] shows that not only is there no 
substantial population in any redshift bin z > 0.6 of low- 
luminosity quasars with masses higher than the rest of the 
quasar population, but also that the highest mass quasars 
at every redshift are restricted to luminosities well above 
the SDSS detection threshold. The luminosity distribution 
of high-mass quasars in four representative redshift bins is 
shown in Figure [7] 

At almost every mass and redshift the peak quasar num- 
ber density occurs at a luminosity above the SDSS detection 
threshold (Figure[T|. In Paper I, we defined the SEB by look- 
ing at the decline from that peak on the high-luminosity end. 
Here we similarly define a boundary on the low-luminosity 
end as the point at which the number density has dropped 
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Table 3. Best-fitting linear HMLLB slopes at 0.2 < z < 2.0. 
There is no HMLLB above the SDSS detection threshold at z < 
0.6. Where an HMLLB exists, its slope seems not to be a constant 
fraction of L Edd , but rather at or below that of the SEB at the 
same redshift. The range of L/L^dd from the SEB to the HMLLB 
at the high-mass end is also included. 



z Line Slope SEB slope L/L^dd range (%) 



0.2-0.4 H/3 none 0.37 ± 0.02 1-8 

0.4-0.6 H/3 none 0.45 ± 0.03 1-11 

0.6-0.8 H/3 0.60 ±0.19 0.60 ± 0.06 1-20 

0.6-0.8 Mgll 0.67 ±0.10 0.61 ±0.10 3-32 

0.8-1.0 Mgll 0.54 ±0.03 0.67 ± 0.09 3-25 

1.0-1.2 Mgll 0.35 ±0.03 0.67 ± 0.05 2-25 

1.2-1.4 Mgll 0.44 ±0.07 0.73 ± 0.05 2-31 

1.4-1.6 Mgll 0.47 ±0.02 0.68 ± 0.08 3-33 

1.6-1.8 Mgll 0.44 ±0.04 0.50 ±0.10 3-25 

1.8-2.0 Mgll 0.31 ±0.03 0.42 ± 0.06 3-24 



to 20% of peak. Note that because the peak number density 
is not constant across mass bins, this boundary is different 
from a 20% contour on the entire locus. Contours of quasar 
number density in the M — L plane also exclude a high-mass, 
low-luminosity region (Figure [T}. 

We consider a high-mass, low-luminosity boundary 
(HMLLB) to be present when for at least the highest 0.5 
dex in mass, the boundary lies significantly above the SDSS 
detection threshold. Best-fitting lines to the HMLLB, where 
present, are summarized in Table [3] Poisson uncertainties on 
the number density of objects at the boundary suggest a rel- 
ative uncertainty between the boundary at different masses. 
We use these relative uncertainties in the determination of a 
best-fitting line. At each sufficiently high redshift (z > 0.6), 
there is a HMLLB present. Below this redshift, the peak 
number density lies near the SDSS detection limit, so it 
is possible that a HMLLB would exist with improved de- 
tection. The HMLLB is not at a constant fraction of L E dd 
(slope of 1.0), but rather at a slope at or below that of the 
sub-Eddington boundary. 

At 0.6 < z < 0.8, where there are both Mgll and H/3 
masses available, the HMLLB have similar slopes but start 
at a ~ 0.2 dex lower mass in the Mgll quasar locus than in 
the H/3-based locus. Onken & Kollmeier (2008) suggest that 
Mgll masses may need to be recalibrated against H/3 masses, 
and Risaliti, Young, & Elvis (2009) compare H/3 and Mgll 
masses to suggest \og[M B H(H (3)] = 1.8xlog[MsH (Mgll)]- 
6.8.. Applying this correction would produce a Mgll-based 
HMLLB slope of 0.37 ±0.06 at 0.6 < z < 0.8. Both this cor- 
rected slope and the uncorrected HMLLB slope of 0.67±0.10 
are consistent with the H/3-based slope of 0.60 ±0.19. While 
the SEB at this redshift implies the Mgll mass correction 
may not be required for the most massive and most lumi- 
nous objects at each redshift, the HMLLB only provides an 
inconclusive test of these corrections at high mass and low 
luminosity. 

Possible explanations for the dearth of high-mass, low- 
luminosity quasars can be divided into three categories: 

(i) The SDSS pipeline might not select such objects as 
quasars, either due to an error in selection or because there 
is a high-mass, low-luminosity boundary (HMLLB) at which 
quasars gradually transition to an object classified differ- 
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Figure 8. Colors of the 15 quasars (red) at 1.4 < z < 1.6 with 
M > 10 9 5 Mq and L < 0.025L Edd compared to all SDSS quasars 
(black). These high-mass, low-luminosity quasars do not all lie 
along a common photometric selection boundary but are comin- 
gled with other SDSS quasars. 

ently (such as a Seyfert galaxy) but with a similar spectral 
energy distribution to a quasar. 

(ii) Their masses might be systematically underestimated 
or luminosities systematically overestimated. 

(iii) There might be a high-mass, low-luminosity bound- 
ary at which quasars very rapidly transition to an object 
with a very different spectral energy distribution outside 
the SDSS locus, spending little time as a transitional object 
similar to a quasar. 

We will consider the first two possibilities below and the 
third in § H 

SDSS selection for well-measured objects above the de- 
tection limit consists of identifying the object from its mag- 
nitudes in the five SDSS spectral bands. If there is an HM- 
LLB at which SDSS mistakenly ceases to select objects as 
quasars, this would imply that quasars at the HMLLB are 
also at a boundary of the SDSS quasar locus. The sharpest 
HMLLB (with the fewest points lying below the HMLLB 
and the best-determined slope) occurs in the 1.4 < z < 1.6 
redshift bin, so it is the best candidate for a selection effect. 
At that redshift, there are 15 quasars with M > 10 9 5 M Q 
and L/LEdd < 0.025. In Figure [8] the colors of these 15 ob- 
jects are shown in four slices of the five-dimensional color 
space used by SDSS in its photometric selection. Their lo- 
cations are not tightly clustered in color space, although 
a larger fraction lie at the outer edges of the SDSS loci 
than would be expected if they were drawn randomly from 
the SDSS quasar catalog. However, the HMLLB does not 
lie along a photometric selection boundary, nor are there 
additional selection criteria that these objects are close to 
URichards et al. 20020 . SDSS select ion is not responsible for 
the HMLLB. 

While a mass-dependent systematic error in luminos- 
ity determination seems plausible, masses at the HMLLB at 
lower redshift are at the low-mass end of the quasar sample 
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Figure 9. Spectrum of SDSS J171737.29+535647.3, a quasar ly- 
ing on the HMLLB at z = 1.46. The FWHM of Mgll as deter- 
mined by Shen et al. (2008) is indicated in red. The FWHM re- 
quired to increase the virial mass estimate by 0.5 dex is indicated 
in blue. 

at higher redshift. As discussed in Paper I, when quasars de- 
termined to have identical masses are at the low-mass end 
of the quasar sample, they reach Lsdd but do not surpass 
it. Therefore, any systematic over-estimation of quasar lumi- 
nosity of the sort required to explain the HMLLB either also 
requires that no high-redshift quasar be near Lsdd or that 
the HMLLB is caused by a coincidence between redshift- 
dependent and mass-dependent systematic errors that al- 
most entirely cancel for quasars at lower masses in every 
redshift bin while producing the HMLLB at higher masses. 
Such fine-tuning seems highly implausible, although it can- 
not be ruled out from direct observational evidence. 

A systematic error in mass determination at the HM- 
LLB requires virial masses to be incorrect for faint quasars 
with very broad emission lines. For Mgll in particular, the 
wings of wide emission lines with low equivalent width might 
be difficult to discern from slightly redder Fell lines. To pro- 
duce the HMLLB, at 1.4 < z < 1.6 high-mass and low- 
luminosity objects would need to have their masses under- 
estimated by ~ 0.5 dex, if the highest-mass objects have 
identical mass at every luminosity (Figure [1} . If turnoff is 
slow and the highest mass objects are more massive at low 
L then a greater factor would be required. Because virial 
masses depend upon the square of the H/3, Mgll, or CIV 
full- width half-maximum (FWHM), an error of ~ 0.5 dex in 
mass would require an error of ~ 0.25 dex in FWHM. 

Figure [9] shows the spectrum of a quasar lying near 
the 1.4 < z < 1.6 HMLLB. This spectrum has the median 
signal-to-noise ratio for the 1.4 < z < 1.6 HMLLB of 3.36, 
while quasars on the HMLLB range in signal-to-noise from 
2.26 to 6.29. The FWHM measured by the SDSS pipeline is 
shown in red. To increase the virial mass estimate by 0.5 dex, 
the FWHM would need to instead be of the size indicated 
by the blue line, which is clearly too broad and is closer to 
the entire width of the line itself than to its FWHM. 

There is an window of redshift at 0.6 < z < 0.8 in 
which both H/3-based and Mgll-based mass estimates are 



available. If the Mgll-based estimates are systematically too 
small near the HMLLB due to Fe contamination, the H/3- 
based estimate for the same objects should be systemati- 
cally larger. For all 3505 objects at 0.6 < z < 0.8, the mean 
Mm 9 — Mh/3 = —0.06 dex, with a standard deviation of 0.28 
dex. Of the 56 objects closest to the HMLLB (as determined 
by Mgll masses), the mean Mug — Mh/3 = 0.05 dex smaller 
than those of Mgll, with a standard deviation of 0.44 dex. 
Because these are at the high-mass tail of the Mgll mass 
distribution, they are more likely to be objects where the 
Mgll mass is larger than the H/3 mass if their statistical un- 
certainties are uncorrelated. So, the atypically large mean 
Mm 9 — Mh/3 for this sample should be expected. However, 
neither the entire population nor only quasars near the HM- 
LLB have a difference close to the the ~ 0.5 dex correction 
required to remove the HMLLB. 

Finally, since virial mass estimation is calibrated from 
just ~ 30 objects (Vcstcrgaard & Peterson 2006), the rela- 
tion may be poorly calibrated for unusual objects, perhaps 
including those near the HMLLB. In particular, it is possible 
that the line FWHM and continuum are measured correctly 
but virial masses are systematically low near the HMLLB. 
This possibility cannot be ruled out directly from available 
data. Virial mass estimates can be very nearly derived the- 
oretically under the two assumptions that (1) peculiar ve- 
locities in the broad-line region are nearly virial and (2) 
the luminosity is nearly L oc R 2 . The successful calibra- 
tion of virial mass estimates against reverberation mapping 
suggests that these two assumptions hold for most quasars. 
Breakdown in virial mass estimation near the HMLLB would 
be an indicator that the properties of quasars change near 
the HMLLB. 

In conclusion, the HMLLB is not caused by (1) erro- 
neous SDSS selection or (2) errors in determining spectral 
parameters used in virial mass estimation. Incorrect lumi- 
nosity determination near the HMLLB would require a fine- 
tuning between redshift-dependent and mass-dependent sys- 
tematic errors. The most likely explanation for the HMLLB 
appears to be a transition in the properties of the quasar, 
which may be manifested either as a rapid change in the 
quasar SEB (leading to an object that no longer resembles 
a quasar in the five SDSS color bands) or as a breakdown in 
virial mass estimation near the HMLLB (e.g., if the broad 
line region is no longer predominantly virial). 

4.2 Line Ratios 

Different emission lines come from transitions with different 
ionization potentials, and are interpreted as being emitted 
at different radii from the central black hole (Peterson 2008). 
In particular, Mgll is emitted furthest from the central black 
hole and CIV nearest. Changes in quasar structure near 
turnoff, either out in the broad-line region (BLR) or closer in 
at the optical/UV emitting part of the accretion disk, may 
show up as unusual line ratios. The redshift ranges where 
we can see either the Mgll and H/3 or Mgll and CIV lines in 
the same spectrum allow us to test this idea. In Figure 10 
we show these two line ratios as a function of mass. 

For 0.6 < z < 0.8, the H/?/Mg/small II ratio has a 
similar distribution at all black hole masses, although the 
very highest-mass objects (blue, Figure [10)) might hint at a 
slight deviation (3% chance to be identical, A Kolmogorov- 
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Figure 10. Distribution in three mass bins of the flux ratio of 
the H/3 to Mgll (top) and Mgll to CIV (bottom) emission lines 
at redshifts where both lines are visible in SDSS spectra. At top 
(z ~ 0.6), the boundaries between bins are log M/Mq =8.0 and 
9.0, while on the right (z ~ 2.0), the mass boundaries are at 
log M/Mq = 9.0 and 9.5. The high-mass distribution is drawn in 
blue while the low-mass distribution is drawn in red in each panel. 
The indicated uncertainties are Poisson errors and are likely un- 
derestimates. 

Smirnov test). This similarity suggests that, high-mass and 
low-mass quasars at z ~ 0.6 have both similar BLRs and 
ionizing continuum shapes. There is no difference in either 
the virial mass or bolometric luminosity estimation tech- 
niques between high-mass and low-mass quasars, and this 
similarity is at least partial confirmation that this might be 
reasonable. 

A majority of the 723 quasars at M > 9.4 will turn off 
by 2 = 0.4, i.e., within their next ~ 1.4Gyr. Since we see 
no difference in the H/3/MgII line ratio, this either means 
that the process by which they cease luminous accretion 
does not substantially change this ratio or that the turnoff is 
sufficiently rapid that we see no signature of it in Figure [lOl 
Since there are ~ 400 quasars that turn off within ~ 1.4 Gyr 
and there is no substantial population in the tails of the line 
ratio distribution, if the turnoff process involves one of these 
two emission lines substantially disappearing, we have seen 
no examples and therefore the turnoff process must require 
no more than tens of Myr to complete. Alternatively, this 
places no limits on the timescale for a low-redshift turnoff 
process involving no changes in these line ratios. 

At 1.8 < z < 2.0, the Mgll to CIV line ratio changes 
with mass, with the Mgll line being stronger at high black 
hole mass. One possibility is that massive black holes tend to 
have strongly asymmetric CIV lines (|Shen et al. 2008(1 . and 
therefore the line flux is underestimated by a Gaussian fit. 
However, the CIV line has a higher ionization potential. As 
a result, perhaps we could interpret Figure [10] as showing 
that massive black holes have a different ionizing continuum 
than less massive black holes near z ~ 2.0. This might give 
us a hint as to the turnoff mechanism for massive black holes 
at z ~ 2.0. 



It is curious that the Mgll to CIV ratio changes with 
mass 1.8 < z < 2.0 while at 0.6 < z < 0.8 the H/3/MgII ratio 
is not affected. Because the three lines come from different 
ionization potentials, this apparent discrepancy might indi- 
cate that the turnoff process acts differently on the broad- 
line region and the ionizing continuum. Alternatively, it 
might indicate that the turnoff mechanism is changing be- 
tween these two redshift ranges. Additional spectroscopy in 
the infrared and ultraviolet is possible and will be able to 
determine which interpretation is correct. 



5 A SYNCHRONIZATION PUZZLE 

In § [3] it was shown that the number density of quasars 
of a given Mbh declines with a mass-dependent e-folding 
time between 0.7 and 3 Gyr for quasars with an Mbh of 
1O 9 M - 10 10 M Q (Table The SDSS catalog includes 
much of the Northern hemisphere and so includes quasars 
almost diametrically opposed in the sky. At a redshift of 
2, they lie in host galaxies that had not been causally 
connected since inflation, while even at much lower red- 
shift it is difficult to believe galaxies a few Gpc apart 
would strongly influence each others' development. Yet, 
quasars with Mbh ~ lO lo Af0 in such galaxies turn off syn- 
chronously to within 700 Myr. 

Could this be resolved if virial mass estimates are poorly 
calibrated at high mass? While the details of quasar turnoff 
do rely upon high-mass objects, there is a similar synchro- 
nization in the sub-Eddington boundary (Paper I). This 
boundary restricts quasars to a luminosity below LEdd at 
most masses and every redshift, with the restriction stronger 
with increasing mass. However, for every mass there is a red- 
shift above which there is a substantial population accret- 
ing near LEdd- This again presents a synchronization prob- 
lem; for example, a substantial population of quasars with 
Mbh ~ 10 s are at Eddington at z ~ 0.6, but none are at 
Eddington at z ~ 0.4, less than 2 Gyr later. So it is not just 
turnoff that is synchronized for quasars at a given Mbh but 
rather most of their luminous accretion. Because the loca- 
tion of the sub-Eddington boundary declines in luminosity 
over time at a given mass, its sharpness is further evidence 
of this universal synchronization of quasar accretion. 

Each quasar lies within a host galaxy, and it is widely 
assumed that the dynamics of quasar accretion are con- 
trolled by the dynamics of its host, since the host galaxy 
appears to be the only fuel source available. Turnoff, then, 
would occur when the galaxy no longer can sufficiently fuel 
its central black hole, a condition presumably related to 
the evolution of the host galaxy. Further evidence for co- 
cvolution of the quasar and the host galaxy comes from the 
black hole mass - stellar velocity dispersion (Mbh — cr) re- 
lation i|Ferrarese fc Merritt 20001 IGebhardt et al. 2000[> as 
well as because quasars of a given mass tend to also lie in 
host galaxies of a common mass ( Fcrrarcsc & Merritt 2000 
IGebhardt et al. 2000|> . 

However, host galaxies of a given mass are not believed 
to be as well synchronized as required by quasar accretion 
or turnoff. The Millenium simulation (Springcl et al. 2005) ) 
continues to produce new galaxies at low redshift, while 
galaxy merger rates are believed to only decline by 
a factor of ~ 10 between a redshift of 2 and to- 
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day QHopkins et al. 2009| l, equivalent to a characteristic e- 
folding timescale of ~ 5 Gyr. In short, if the dynamics of 
quasars are only as synchronized as those of their host galax- 
ies, our current theoretical understanding of galactic dynam- 
ics is incompatible with the data presented here. 



6 DISCUSSION 

Here and in Paper I , we have considered what the quasar lo- 
cus in the mass-luminosity plane can teach us about quasar 
accretion and turnoff. Quasars have been modeled as very 
simple objects, in part because there has not yet been suf- 
ficient data to require more complex assumptions. In Paper 
I, we falsified the assumption that quasars at all combina- 
tion of mass and redshift were near their Eddington lumi- 
nosity while accreting luminously. With a combination of 
the large quasar catalog from the Sloan Digital Sky Survey 
(SDSS) and virial mass estimation, we can now evaluate 
several other common assumptions. 

It is widely assumed that the dynamics of quasar turnoff 
are closely linked to the dynamics of the host galaxy. How- 
ever as shown in § O accretion rate and turnoff are highly 
synchronized amongst quasars of a common mass. Either 
galaxies must be much more highly synchronized than cur- 
rently believed or this synchronization must have an alter- 
nate source. Galaxies continue to merge and virialize at later 
times than their quasars are predominantly turning off, and 
the decline in galaxy merger rates is believed to be much 
more gradual than that observed for quasars. So, why are 
quasars so well synchronized? One possibility is that the 
seeding mechanism comes not from the dynamics of their 
host galaxy but rather from some mechanism with a cos- 
mic clock. The ~ 700 Myr timescale for quasar turnoff at 
Mbh ~ 1O 1O M0 corresponds to the Hubble time at z ~ 8, 
and extrapolating the slope of the e-folding time as a func- 
tion of mass (§ [3| to the highest-mass quasars would yield 
the Hubble time at z ~ 20. It is plausible that the seeds of 
these massive quasars turn on during the epoch of reion- 
ization and, if the seeding mechanism is one that yields 
a nearly-fixed quasar life span, turnoff would be synchro- 
nized to roughly the observed precision. This would require 
a seeding mechanism only weakly linked to the dynamics 
surrounding galaxy, but rather triggered by cosmic dynam- 
ics. 

It is also widely assumed that most quasars are near 
Eddington while accreting luminously. The presence of a 
sub-Eddington boundary at which high-mass quasars can- 
not reach their full Eddington luminosity directly contra- 
dicts this claim and slows the maximum allowed growth 
rate. However the boundaries of the quasar locus are very 
sharp, and quasars of a given mass are highly synchronized 
in both the maximum luminosity they can attain and in their 
turnoff time. In the absence of techniques for mass estima- 
tion of high-redshift quasars, it was widely assumed that a 
segregation of quasars by mass was equivalent to a segre- 
gation by luminosity. Since individual quasars are variable, 
this cannot be entirely true, but there is certainly a corre- 
lation between quasar mass and luminosity as evidenced by 
the quasar distribution in the mass-luminosity plane (Figure 
[TJ. Further, the high- mass, low- luminosity boundary sug- 
gests that even if high-mass quasars cannot reach LEdd, 



luminous accretion is rare below a few percent of Edding- 
ton as well. This also implies that the large Milliquas sam- 
ple l|Richards et al. 2009[) of ~ 10 6 new SDSS quasars se- 
lected photometrically but not bright enough for SDSS spec- 
troscopy consists predominantly of quasars at lower Mbh 
than those in the SDSS spectroscopic quasar catalog rather 
than a mixture of quasars at low Mbh and quasars at higher 
Mbh and very low Eddington ratio. 

Finally, it is widely believed that a typical SMBH 
will have many periods of luminous accretion separated 
by dormant periods. In part, this is because the Salpeter 
time for black hole accretion at LEdd. is only ~ 10 8 yr 
( |Natarajan fe P ringic 1998]) , and a quasar could not grow at 
LEdd for several Gyr. However, the sub-Eddington boundary 
restricts more massive quasars at every redshift to 0.2 Lsdd 
or below (with the peak number density below 0.1 LEdd), 
raising the possibility that rather than accreting luminously 
for less than 1 Gyr, quasars might accrete luminously, but 
sub-Eddington, in one long, continuous burst for several Gyr 
from when they first turn on until turnoff. Certainly the 
sub-Eddington boundary requires that quasars accrete lu- 
minously far longer than previously believed. 

When quasars of a given mass no longer reach a given 
Eddington ratio, the sub-Eddington boundary never moves 
back to a higher luminosity. When quasars of a given mass 
begin to decline in number density, their turnoff is synchro- 
nized and sharp, and the population does not reappear at 
lower redshift. While this is insufficient evidence to con- 
clude that turnoff is a permanent process for any individual 
quasar, a quasar lifecycle consisting of one long period of 
accretion followed by permanent turnoff is consistent with 
all of these observations. 

As shown in § |4j the sharpness of the high-mass, low- 
luminosity boundary (HMLLB) appears to require either 
that virial mass estimation is wrong near the HMLLB or 
that quasars at the HMLLB undergo a rapid transition to 
an object with a very different spectral energy distribu- 
tion. One possibility would be transition to RIAF accretion 
Parayan fc McClintock 2008) . However, the HMLLB does 
not occur at a constant Eddington ratio as accretion disk 
theory would suggest, but rather has a sub-linear depen- 
dence on mass. The HMLLB is also evolving with redshift: 
the HMLLB for high-mass quasars is at a higher redshift 
than lower-mass quasars. 

In this paper, we have shown several new puzzles pre- 
sented by the quasar mass-luminosity plane in addition to 
the sub-Eddington boundary discussed in Paper I. These 
new puzzles which cannot be explained by selection effects 
are further evidence that supermassive black hole accre- 
tion and turnoff is more complex than previously believed. 
Rather, many of the common assumptions made in quasar 
toy models are poor matches to the quasar locus in the M—L 
plane. These new boundaries indicate both that a more nu- 
anced quasar model is required and that the quasar mass- 
luminosity plane now contains enough detail that it might 
be able to provide strong tests for such new models. 
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